Quantum noise memory effect of multiple scattered light 
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We investigate frequency correlations in multiple scattered light that are present in the quantum 
fluctuations. The memory effect for quantum and classical noise is compared, and found to have 
markedly different frequency scaling, which was confirmed in a recent experiment. Furthermore, 
novel mesoscopic correlations are predicted that depend on the photon statistics of the incoming 
light. 
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Light propagating through a disordered distribution of 
strongly scattering particles induces an apparently ran- 
dom intensity pattern known as a volume speckle pat- 
tern. Despite the apparent randomness, such intensity 
speckle patterns can posses strong temporal and spatial 
correlations, giving rise to the so-called memory effect "l^. 
The memory effect manifests itself as strong correlation 
between different spatial directions or frequency compo- 
nents of the speckle pattern. After its demonstration in 
early experiments ^ ^ , the memory effect has been em- 
ployed as a sensitive technique of measuring the diffusion 
constant of light propagation For extremely strong 
scattering, close to the transition to Anderson localiza- 
tion, classical intensity correlations can be strongly en- 
hanced, which are referred to as mesoscopic correlations 
0, fij . Recently, it was shown that speckle correlations 
also exist when recordingthe fluctuations of light instead 
of merely the intensity 7]. Furthermore, strong spatial 
quantum correlations were predicted depending on the 
quantum state of light incident on the medium 8]. In 
the present Letter we derive the noise correlation func- 
tion, which was measured in 7] , for an arbitrary quantum 
state of light. Our work demonstrates that the quantum 
fluctuations of a volume speckle pattern give rise to a 
strong memory effect that behaves different than classical 
fluctuations and also different than the classical memory 
effect for intensities. 

We study the propagation of photons through a mul- 
tiple scattering random medium. Using the quantum 
model for multiple scattering we discretize the sys- 
tem into N different input modes and output modes. Let 
h^'' denote the photon number operator for light at the 
frequency lu propagating from input channel a to out- 
put channel b. The fluctuations are quantified by the 



where the quantum mechanical expectation value will be 
evaluated for different quantum states of light. Further 
details about this model are given in The quantum 
noise frequency correlation function for light propagating 
from input channel a to output channel b is defined by 
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This function measures the correlation between the fiuc- 
tuations in a specific output direction for a frequency off- 
set of Auj, and is the generalization of intensity correla- 
tion functions to the case of photon number fluctuations. 
The bars above the variances denote averaging over en- 
sembles of disorder. We consider an arbitrary quantum 
state of light coupled through channel a characterized by 
an averaged number of photons (nJJ,) and a Fano factor 
Fa — (An^)^/ (n°}. The photon number fluctuations of 
light transmitted to channel b is given by Q 

[i^<'f]QN = T^;' + iK) (.Fa - l)iTff, (2) 

where T^'' is the intensity transmission coefficient from 
channel a to channel b. Here the subscript QN refers to 
quantum noise of the relevant quantum state. In the sit- 
uation where classical noise (CN) is dominating, which, 
e.g., is the case when modulating a laser beam, the vari- 
ance of the fluctuations is proportional to the squared 
average number of photons |3 

[iAnZ'>r]^^^{KfiTff. (3) 

Consequently, the correlation functions for quantum and 
classical noise are given by 
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where we have assumed {fiuj) = {n^+Aui) corresponding 
to the situation where the input number of photons is 
kept constant when varying the optical frequency, which 
was the case in the experiment in |7j . Note that we have 
skipped the indices ab on the transmission coefficients in 
order to simphfy notation. 

In the lowest order approximation, the electric field 
transmission coefficient ti^ can be described by a circular 
Gaussian process, where T,,. = \tujf. For such Gaussian 
random variables, we have |inl | 

tl . . . i^ife+l ■ ■ ■t2k— ^ titp X ^2^q X ... X t^tr, (5) 
7r 

where the summation is over all k\ permutations of the 
indices. Based on this theorem, we can evaluate the mo- 
ments of that are relevant for Eqs. (|4a() and 14bl) : 
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T.^ = nlT, 
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First, for the quantum noise, let us restrict to the special 
case of shot noise (SN), which corresponds to = 1. 
Hence, the correlation functions in Eqs. (|4a() and (|4b|) 
are given by 

2 

= /(Ac.), (7a) 
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is the function that describes the frequency decay of the 
intensity correlations with ojd = D/2L^. 

The correlation functions for shot noise and technical 
noise are plotted in Fig. ^ as a function of frequency 
offset Aw. We observe that the noise correlations are 
much stronger in the case of classical noise relative to 
shot noise, thus the correlation at Aw = is a factor 5 
higher. Thispronounced difference was verified experi- 
mentally in [Tf . We note that the correlation function for 
shot noise is identical to what one would get if recording 
the intensity, the latter being independent of the actual 
quantum state. Notably, the quantum fluctuations pro- 
vide an independent measure of frequency speckle corre- 
lations, and as will be seen in the following, the noise cor- 
relations depend strongly on the quantum state of light 
used in the experiment. 



FIG. 1: Noise correlation function for shot noise (full curve) 
and classical noise (dashed curve) as a function of frequency 
offset. 



In the general case of an arbitrary quantum state char- 
acterized by the Fano factor Fa , the correlation function 
is given by Eq. H4a(l . This expression contains terms of 
various orders in Ti^, where in general Ti^ <C 1, i.e. Eq. 
(|la|l can be expanded in orders of the intensity trans- 
mission coefficient. In this expansion of the products of 
transmission coefficients also higher-order contributions 
must be included that originate from deviations from 
Gaussian statistics of the transmission coefficients. The 
lowest order contribution of such mesoscopic correlations 
provides a correction term in the product of two trans- 
mission coefficients 

2 2 3 

TujT^+Aiu = T^ + \tZtiu+Auj\ + "^^5'(Aw)r„ , (9) 

where L is the thickness of the random medium, £ is 
the transport mean free path, and we have defined the 
function 



UJD smh{ J Auj/lud) ~ sin{ J Auj / lu d) 
g{Auj) = — X ^ ^ i-. 
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Expanding Eq. H4a|) to first order leads to 

C^.'^iAcj) « CUAu;) + CiiiAu;), (If) 

where 
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The dominating term in the expansion (C^^^) is simply 
equal to the contribution obtained with shot noise. De- 
viations from shot noise behavior is observed when using 
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Ti^ as a function of frequency offset for three different quan- 
tum states of light: Fock state Fa = (sohd line), coherent 
state Fa — 1 (dashed line) and thermal state Fa — 2 (dotted 
line). The ratios of mean free path to the sample thickness 
are i/L = 1/3 (bold lines), £/L = 1/4 (medium lines), and 
£/L = l/5 (thin lines). 



quantum states of light different from the coherent state, 
i.e. Fa ^ 1, which gives rise to the second term in Eq. 
(|12b|l . The quantum correction competes with classical 
mesoscopic correlations (first term in Eq. (|12b|l ). which is 
in contrast to the dominating quantum corrections found 
in the fluctuations of the total transmission and reflection 
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In Fig. 12 we plot the second-order noise correlation 
function C^^(Acj) for three different quantum states of 
light, corresponding to = 0, Fa = 1, and Fa = 2. 
These Fano factors can be achieved with single-mode 
Fock states, coherent states and thermal states, respec- 
tively Q . Figure [21 also indicates the correlation func- 
tion for different ratios £/L, where £/L ^ 1 for optical 
media where multiple scattering dominates. The classi- 
cal mesoscopic correlations, which also can be extracted 
from intensity measurements, are obtained for Fa = 1. 
In the limit Aoj/ujd — > 0, the classical correlation func- 
tion diverges, which is a consequence of the plane- wave 
approximation and is suppressed when including flnite- 
width beams .12]. This sensitivity to the width of the 
beam only plays an important role for the correlations at 



Aui/lod ~ 1. We observe from Fig. I^jthat either posi- 
tive or negative quantum noise correlations are obtained 
using either super-Poissonian {Fa > 1) or sub-Poissonian 
photons {Fa < 1), respectively. Consequently, the fluctu- 
ations are found to possess novel correlations depending 
on the quantum state of light, which is markedly differ- 
ent from intensity correlations that are independent of 
the quantum state. 

A novel noise correlation function for multiple scat- 
tered light was introduced and evaluated for both classi- 
cal noise and for arbitrary single-mode quantum states. 
Pronounced different correlations were found when com- 
paring classical noise to quantum noise. Including higher- 
order correction terms in an expansion in the transmis- 
sion coefflcient, quantum corrections to the noise corre- 
lation function were predicted that have no analogy in 
classical intensity measurements. 
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